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ABSTRACT 

The  momentum  theory  of  moderately  loaded  propellers  has  been  generalized 
to  account  for  the  effects  of  shear  in  the  oncoming  stream.  A condition  to 
determine  the  distribution  of  thrust  for  optimum  efficiency  is  derived  which 
generalizes  a condition  derived  by  Betz  for  a propeller  in  a uniform  flow. 

This  condition  leads  to  a nonlinear  integral  equation,  and  a computer  program 
has  been  developed  for  solving  It.  Results  for  a propeller  of  a given  radius 
in  particular  wakes  at  different  Reynolds  numbers  are  presented  and  discussed. 
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Dirac  delta  function;  also  symbol  for  first  variation 
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INTRODUCTION 

The  theory  of  propeller  design  is  a we  1 1 -developed  art  that  has  pro- 
gressed for  over  a century  from  axial  momentum  theory  to  general  momentum 
theory  to  blade  element  theory  to  lifting-line  theory  to  lifting-surface 
theory.  It  has  thus  become  possible  to  predict  the  mean  radial  pressure 
distribution  with  various  degrees  of  refinement  for  a given  propeller  as 
well  as  the  harmonic  content,  and  also  to  determine  the  mean  load  distribu- 
tion for  optimum  efficiency. 

Despite  the  sophisticated  mathematical  treatments  presented  by  many 
authors,  all  theories  thus  far  developed  have  been  based  on  the  existence 
of  a velocity  potential.  Propellers  which  are  placed  near  the  stern  of  a 
ship  or  submarine  operate  in  the  wake  of  the  vehicle  where  the  oncoming 
stream  contains  vorticity  and  this  is  usually  taken  into  account  by  para- 
metrically varying  the  wake  velocity  over  the  plane  of  the  propeller  and 
retaining  the  assumption  of  potential  flow.  This  procedure  can  be  justified 

only  if  the  vorticity  is  weak  (see  Sears,'  p.^2).  Thus,  the  lifting  surface 

2 3 k 

theories  presented  by  Morgan  and  Wrench,  Yamazaki,  and  Greenberg,  are  all 
based  on  this  assumption,  and  Greenberg  goes  so  far  as  to  say  that  unless 
the  assumption  is  made  "..the  problem  becomes  hopelessly  unmanageable."  This 
is  very  likely  to  be  so  for  general  unsteady  lifting  surface  theory.  However, 
it  leaves  unanswered  the  question  as  to  the  suitability  of  the  assumption 
for  an  actual  ship  wake  whose  vorticity  may  or  may  not  be  sufficiently  small 
to  Justify  the  assumption. 

In  those  cases  where  the  vorticity  is  not  weak,  a velocity  potential 
cannot  be  assumed  and  the  presence  of  shear  in  the  oncoming  stream  will  af- 
fect the  flow  field.  Theories  which  account  for  the  shear  have  been  considered 
by  several  authors  for  the  floiv  about  wings,  the  first  such  paper  being  that 
of  Karman  and  Tsien,^  but  the  parallel  development  for  propellers  seems  to 
have  been  ignored  probably  because  of  the  unmanageability  of  the  problem. 


Superior  numbers  in  text  matter  refer  to  similarly  numbered  references 
listed  at  the  end  of  this  report. 
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Nevertheless,  it  would  be  useful  to  have  some  simple  solutions  which 
do  account  for  shear  so  that  comparisons  can  be  made  with  the  weak  shear 
approximation  thereby  allowing  an  assessment  of  the  weak  shear  assumption. 

Two  possibie  "simple"  theories  appear  to  be  possible  in  which  the  effect  of 
shear  is  included. 

The  first  such  is  momentum  theory  which  is  basically  a one-dimensional 
theory  and  for  which  the  wake  profile  can  be  allowed  to  vary  radially  in  an 
arb i trary  way . 

The  second  such  theory  is  a steady  state  lifting-line  theory  in  which 
shear  is  accounted  for  but  the  wake  profile  takes  on  a special  shape  for 
which  the  field  equation  becomes  no  more  complicated  than  that  for  potential 
flow. 

The  analysis  presented  below  is  confined  to  the  first  of  these  simpli- 
fied theories,  but  it  is  intended  to  develop  the  second  in  a subsequent  re- 
port. In  Section  I the  equations  for  the  axial  momentum  theory  are  developed. 
In  Section  II  the  equations  for  the  general  momentum  theory  are  developed  and 
the  condition  of  Betz  for  the  optimum  distribution  of  thrust  is  generalized 
to  an  arbitrary  wake  distribution.  In  Section  III  some  numerical  evaluations 
are  presented  for  a measured  wake  and  comparisons  made  with  results  obtained 
simply  by  varying  the  wake  parametrically  (weak  shear  approximation). 

The  analysis  will  generally  follow  the  article  by  Glauert^  for  the 
momentum  theory  of  a propeller  in  a uniform  stream.  This  article  synthesizes 
the  main  developments  of  the  theory  which  have  changed  very  little  since  it 
was  written  (but  see  Reference  7 for  a recent  development). 

This  research  was  sponsored  by  the  David  Taylor  Naval  Ship  Research 
and  Development  Center,  General  Hydromechanics  Research  Program,  under  Con- 
tract N00014-77-C-0297,  DL  Project  4512/023. 
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I.  THE  AXIAL  MOMENTUM  THEORY 


We  consider  the  Euler  equations  of  motion  for  an  incompressible 
frictionless  fluid.  At  x=0  , there  is  a propeller  to  be  modelled  by  an 
actuator  disc  that  can  sustain  a pressure  jump  Ap=  p(x=0+) “p(x=0-) . The  axial 
momentum  equation  in  cylindrical  coordinates  (x,r)  for  axially  symmetric  flow 
i s 

du  . du  . 1 dp  1 

p 5x  p 


dx 


+ + = i Ap(r)6(x)  ; Ap(r)  = 0 for  r > 


(1.1) 


Here 


u 

V 

P 

P 

6 


is  the  axial  component  of  velocity 
is  the  radial  component  of  velocity 
is  the  pressure 
is  the  fluid  density 
is  the  Dirac  delta  function 


The  equation  of  continuity  is 


(ru)  f Ip  (rv)  . 0 


(1.2) 


Now  assume  there  exists  an  oncoming  stream  U(r)  , the  so-called  nominal 
wake,  which  has  been  generated  by  a body  somewhat  ahead  of  the  actuator  disc. 
Since  U(r)  is  taken  to  be  independent  of  the  azimuthal  angle,  we  will  be 
considering  only  the  zeroth  harmonic  of  the  wake  and  the  mean  forces  on  the 
propeller.  The  perturbation  velocity  components  imparted  to  the  fluid  by  the 
u'  and  V'  . Then,  upon  setting 


actuator  disc  are  u'  and  v' 
u = U(r)  + u'(x,r) 

V = v'(x,r) 

and  substituting  into  Eq.(1.2),  we  obtain 


(1.3) 


It  (ru')  + It  (rV)  = 0 


d^  ^ dr 

On  the  other  hand,  upon  substituting  into  Eq.(l.l),  we  obtain 


(1.4) 


■ I du ' , du ' I 1 ^ 1 A / \ 

~ 4p(r)6(x) 


p dx 


(1.5) 
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where  the  term  v'  has  been  ignored  as  being  of  higher  order.  The 

third  term  on  the  ieft-hand  side  of  this  equation  represents  the  effect  of 
the  shear  of  the  oncoming  stream. 

The  continuity  equation  can  be  satisfied  identically  by  defining  a 
Stokes  stream  function  such  that 

rv'  = ; ru'  = - |i  (1.6) 

in  which  case,  Eq.(1.5)  becomes 


du ' , du ' IdU^  1^ 


(1.7) 


At  this  point  we  will  drop  the  primes  on  u'  and  v'  ; since  we  will  rarely 
be  required  to  refer  to  the  original  u and  v , this  will  cause  little 
confusion. 


In  order  to  solve  for  the  velocity  components  everywhere  in  the  field, 
it  is,  of  course,  also  necessary  to  consider  the  radial  momentum  equation 
and  solve  it  simultaneously  with  the  axial  momentum  equation.  However,  for 
our  purposes  we  will  be  content  to  consider  the  overall  momentum  of  the  sys- 
tem and  with  this  view,  only  Eq.(1.7)  will  be  required.  Far  upstream  (x=-«) , 
all  velocity  components  vanish  and  the  pressure  will  be  taken  to  be  zero. 
Thus,  upon  integrating  Eq.(1.7)  with  respect  to  x from  x=-ob  to  x=0-  , 
it  is  found  that 

Uu(O-)  + ^(0-)  + i ^ M0-)+  j p(O-)  =0  . (1.8) 


t 


On  the  other  hand,  if  we  integrate  Eq.(1.7)  from  x=0+  to  x=oo  , and  assume 
that  any  difference  in  pressure  in  the  wake  from  its  value  far  upstream  is 
of  higher  order,  we  obtain 


u[u^-u(0+)]+i[uf  - u^0+)]  + - ^0+)]  = (1.9) 

Upon  forming  the  quantity  ^ by  adding  Eq.(1.8)  and  Eq.(1.9),  it  is  found 
that 

^ = U Ujo  + ~ ^ 'I'co  + U[u(0-)  - u(0+)  ] 

+ i[u^(0-)- u=(0+)]  + i^CHo-)  - iK0+)]  . (1.10) 


*indicated  by  von  Mises  [8]  to  be  of  third  order. 
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Finally,  if  we  integrate  Eq.(1.7)  across  the  disc  from  x=0-  to  x=0+  , 
we  see  that 

u[u(o-)-  u(o+)]+ i[u^(o-)  - u^(o+)]  + -Jr  ^[no-)-  i(0+)]  = 0 (l.ll) 


Subtracting  Eq.vl*ll)  from  Eq.(l.IO),  we  obtain  for  the  pressure  jump 

4a.  u. 


(1.12) 


from  which  the  thrust  of  the  disc  becomes 


T = J ^pdS  = pj  [uu„+^  + i ^ tJdS 


(1.13) 


where  dS  is  the  eiement  of  area  and  the  integration  is  to  be  carried  out 
over  the  area  of  the  actuator  disc. 

This  formula  is  identicai  to  the  resuit  for  the  thrust  as  obtained  from 
ordinary  momentum  theory  with  the  additional  term  proportional  to  ^ that 
accounts  for  the  effect  of  shear. 


At  this  point  we  may  obtain  a relationship  between  the  stream  function 

and  axial  perturbation  veiocity  at  the  disc  and  in  the  wake  to  the  order  of 

linear  theory.  Now  in  linear  theory,  the  pressure  must  descend  to  - p 

1 2 

Just  in  front  of  the  disc  and  rise  to  + jAp  just  behind  the  disc,  so  that 
the  quantity  p(0+)+(p-)  must  vanish.  Thus,  upon  taking  the  difference 
between  Eq.(1.8)  and  Eq.(1.9),  and  ignoring  the  quadratic  terms  in  u^  , we 
obtain 


u + 1.  ^ ^ = ifu  + 1-  M i 1 

"o  rU  dr  ’^o  2L  » rU  dr 


(1.14) 


where 


u(0+)  + u(O-) 


(1.15) 


Upon  eiiminating  u and  u in  favor  of  ik  and  ili  by  virtue  of  the 

O OB  ’o  « 

definition  (1.6),  it  is  easy  to  show  that  Eq.(l.l4)  becomes 


dU  " U 2\dU  " U ’'oo/ 


(1.16) 


This  is  a differential  equation  for  when  i|(^  is  given,  and  its  general 
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solution  is 

where  K is  a constant  of  integration.  However,  for  a propeller  with  no 

boss,  the  axis  r=0  must  be  a streamline  and  hence  both  and  vanish 

O 06 

there;  whence  K=0  . Thus  to  first  order  the  perturbation  stream  function  in 
the  plane  of  the  propeller  is  half  the  stream  function  in  the  ultimate  wake, 
and  hence  the  perturbation  velocities  are  in  the  same  ratio.  This  relation- 
ship is,  of  course,  the  same  as  had  previously  been  deduced  for  a propeller 
in  a uniform  stream. 

Now  reconsider  Eq.(l.ll)  and  let 

6u  = u(O-)  - u(0+) 

Lif  = iKO-)  - '|r(0+) 

Upon  ignoring  the  quadratic  terms  and  utilizing  the  definition  (1.6),  it 
is  easy  to  show  that  this  becomes 

- U ^ (M)  + ^i|i  = 0 . (1.18) 

This  is  a differential  equation  for  the  jump  in  the  stream  function  across 
the  disc.  Its  general  solution  is 

All  = KjU  . (1.19) 

where  Kj  is  a constant  of  integration.  For  a propeller  with  no  boss,  Aili  must 
vanish  at  r=0,  while  for  a propeller  with  a boss,  A\|i  must  vanish  at  the  radius 
of  the  boss.  In  either  case,  Ki  must  vanish,  and  hence  we  have  proven  that  to 
first  order  there  is  no  jump  in  the  stream  function  across  the  actuator  disc, 
that  is,  no  mass  fiow  is  created  at  the  disc,  and  as  a consequence  of  the 
definition  (1.6),  the  axial  perturbation  velocity  is  also  continuous  across 
the  disc,  a result  that  is  the  same  as  that  for  a propeller  in  a uniform  flow. 

Finally,  we  may  cast  Eq.(1.13)  for  the  thrust  into  another  form  by  con- 
sidering that  dS=2nrdr  and  integrating  the  shear  term  by  parts.  Then,  by 
virtue  of  the  definition  (1.6),  the  thrust  becomes 

= P , (1.20) 

where 

\ 
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V = 2U-  Uj^  (1.21) 

and  Uq  is  the  velocity  of  the  wake  at  the  outer  radius  R of  the  propeller 
disc.  Equation  (1.20)  is  identical  in  form  to  the  expression  for  the  thrust 
without  shear  except  that  the  oncoming  stream  velocity  is  replaced  by  V , 
a quantity  that  will  become  negative  whenever  U(r)<  U|^/2  . 

It  is  necessary  to  point  out  that  the  expression  (1.20)  is  the  thrust 
developed  on  the  propeller  because  it  is  operating  in  a wake  of  velocity 
U(r).  In  other  words,  it  is  the  thrust  that  develops  on  a propeller  tested 
behind  a screen  wake  in  a water  tunnel  test  facility  (neglecting  wall  ef- 
fects). It  is  not  the  thrust  that  develops  on  a propeller  mounted  behind  a 
ship  or  submarine  with  this  wake  because  the  induced  velocity  and  pressure 
fields  created  by  the  propeller  will  increase  the  resistance  on  the  body  over 
and  above  that  which  would  be  present  if  the  propeller  were  absent.  Thus, 
when  it  is  mounted  on  the  ship  or  submarine,  the  propeller  must  not  only 
overcome  the  body  resistance  but  also  the  additional  resistance  induced  on 
the  body  by  the  propeller.  Hence,  Eq.(1.20)  represents  only  the  apparent 
thrust  of  the  propeller  and  for  this  reason  it  is  called  Tg  . Determination 
of  the  induced  resistance  on  the  body  is  an  important  element  in  analyzing 
the  propeller  ship  combination,  but  it  is  not  the  problem  that  concerns  us 
here. 

We  now  consider  the  power  required  to  produce  the  thrust  T which 
will  be  called  Pg  . It  is  clear  that  the  power  at  any  annulus  will  be  equal 
to  the  element  of  thrust  times  the  velocity  (U+Uq)  at  which  it  is  being  trans- 
ported summed  over  all  annular  elements.  In  other  words 


^ ; (u+uo)dT 


(1.22) 


Since  to  first  order  u * ^ , then  upon  substituting  for  the  thrust 

o z • 


density  from  Eq.(1.12),  there  is  obtained 


(1.23) 


Finally,  the  apparent  efficiency  is  given  by  the  ratio  of  the  useful 
work  to  the  power  absorbed; 


U T, 

« a 


(1.24) 
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where  is  the  velocity  of  the  propeller  body  combination  (which  unfortunately 
is  not  known  for  the  present  model  since  it  depends  on  the  connection  between 
the  body  wake  and  the  speed  of  advance).  The  appaient  efficiency  can,  paradox- 
ically, be  greater  than  unity  since  we  have  not  debited  the  propeller  with  the 
induced  resistance  on  the  body. 

We  will  now  determine  the  condition  to  be  satisfied  for  the  propeller  of 
maximum  efficiency.  Thus,  for  a given  thrust,  we  wish  to  minimize  the  power 
used  to  generate  it.  In  other  words,  we  wish  to  set  the  first  variation  of 
Pg+X.U^T^  equal  to  zero,  where  \ is  a Lagrange  multiplier.  In  Eq.(l.23),  we  set 

r 

• (1-25) 

o 

Then, 

6T^  = 

3 40  d 

Rf  , 6u^  . r 

2TTp  M U^6u  + 2Uu  6u  + 77  u ^ 6u  - ■= — -r-  f r,  u (r  ) dr 

^JL  cd  08  0B  hoo  ea  2rdr>J  Iso'l.'  1 
o o 

-(u  + f)  r ^ 

o 

+ \U  (Vdu  + u 6u  )}rdr  = 0 . (1.26) 

00  CO  40  eo  ^ 

Consider  the  term 

-2TTp  J (u  + ^ dr  J ri  6u^(ri  )dri 

o o 

Upon  interchanging  the  order  of  integration,  this  becomes 

-2n  f rdujr  / (u(r^ ) + -f)  ^ dr,  . 
or  ^ 

Thus  the  variation  6u^  may  now  be  factored  out  of  all  the  terms  in  Eq. 

(1.26),  and  upon  applying  the  fundamental  lemma  of  the  calculus  of  varia- 
tions and  letting 

u„  = 2aU„  (1.27) 

00  00 

where  a is  the  axial  interference  factor,  we  find  the  condition  for 
optimum  efficiency: 
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3/0  I fU  \ ,3  1 r^  , r, 

-^r)  * Hu-J  ~77r  “d?7R  J ^ ^ V 


l/Ulf.  f 

r^R 


(1.28) 


In  contrast  to  the  case  of  parametric  dependence  on  U , the  axial  inter- 
ference factor,  when  shear  is  included,  depends  on  conditions  at  all  radial 
stations  across  the  actuator  disc  because  of  the  presence  of  the  integral 
terms.  One  way  to  solve  this  integral  equation  for  a is  by  iteration.  The 
first  step  is  to  specify  a value  of  X . Then  ignore  the  terms  involving 
d(U/Uj 

^ ^ . The  resulting  quadratic  can  be  solved  exactly  at  each  radial  sta- 
tion. The  resulting  values  of  a(r/R)  can  then  be  used  in  the  two  integral 
terms  and  the  quadratic  solved  again.  This  procedure  is  repeated  to  con- 
vergence. 

It  will  be  noted  that  it  is  necessary  to  specify  the  Lagrange  multi- 
plier \ , a quantity  that  is  not,  of  course,  known.  The  apparent  thrust, 
on  the  other  hand,  is  presumed  given.  Thus,  once  a(r)  has  been  found  for 
a given  k , the  apparent  thrust  may  be  calculated  from  Eq.(l.20),  the 
apparent  power  from  Eq.(1.23),  and  the  apparent  efficiency  from  Eq.(1.24). 
Hence,  by  varying  \ over  a range  of  values,  a plot  can  be  prepared  of 
T,  vs.  T).  for  optimum  apparent  efficiency  of  a propeller  of  radius  R in 

d d 

a given  wake.  Note  that  these  curves  will  change  if  the  propeller  radius 
changes  so  that  there  may  exist  an  optimum  propeller  radius  especially  for 
wakes  that  are  not  monotonic. 


The  distribution  of  the  axial  interference  factor  as  obtained  by  this 
method  is  not  likely  to  be  realistic  since  rotational  velocities  are  not 
accounted  for,  but  the  curves  of  ideal  apparent  efficiency  vs.  thrust  and 
power  will  be  very  close  to  the  curves  when  rotation  is  accounted  for  since 
it  is  well-known  from  ordinary  momentum  theory  that  rotation  affects  these 
curves  only  slightly. 

The  equations  for  thrust  and  power  may  be  put  into  coefficient  form 
by  defining  a thrust  and  torque  coefficient; 
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K 

« npn"d" 

Here 

n is  the  number  of  revolutions  per  second  = Cl/2n 
d is  the  diameter  = 2R 

We  define  Ja=U  /nd  , then  Eq.(1.20)  becomes  in  coefficient  form 


Equation  (1.23)  becomes 


(1.29) 


(1.30) 


Kn  1 I It  r 1 r/R  _ 

L^(u“  77K  dTTF  I r 


R j R 


(1.31) 


and  the  apparent  efficiency  becomes 

, (K/V) 

^ (VV) 


(1.32) 


Finally  it  should  be  noted  that  the  design  problem  of  finding  u^  , when 
the  load  distribution  ^p(r)  is  given,  can  be  solved  by  casting  Eq.(1.12) 
into  the  form 


P 


Uu. 


X dU  p 
r dr  J 


(1.33) 


which  follows  from  the  definition  (1.6).  This  integral  equation  for  u 

09 

can  be  solved  by  Iteration  in  the  following  way:  Ignore  the  shear  term 
and  solve  the  resulting  quadratic  for  u^  at  each  r . The  integral  can 
then  be  evaluated  and  the  quadratic  solved  again  including  it.  This 
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II.  THE  GENERAL  MOMENTUM  THEORY 

In  the  axial  momentum  theory  as  presented  in  Section  I,  all  rotational 
motions  have  been  ignored.  In  the  general  momentum  theory  we  consider  them, 
and  hence  we  consider  the  radial  momentum  equation,  viz.. 


r dv  w^  dvi  1 ^ 
L'^  dr  " r ^ ^ ~ ~ p ^ 


(2.1) 


Here  w is  the  azimuthai  component  of  velocity.  Now  it  is  known  from  the 
constancy  of  angular  momentum  that  ahead  of  the  propeller  there  is  no  rota- 
tion, while  behind  it  the  propeller  imparts  a rotational  velocity  w to 
the  fluid.  Since  v and  u are  continuous  across  the  propel ier  disc,  we 
obtain  for  the  jump  in  pressure 


J.  dAp  _ 
p dr  r 


(2.2) 


In  a coordinate  system  rotating  with  the  propeller  blades,  the  flow  will  be 
steady.  Hence,  let 


w = -Qr  + w' 


Equation  (2.2)  then  becomes 


(2.3) 


= 2nw'  - 


(2.4) 


Now  let  w'=our  , and  integrate  with  respect  to  r . Then,  if  u)  is  a 

'ft 

slowly  varying  function  of  r , it  is  easy  to  show  that  approximately 


Ap  “ puJ  r=(n  - j) 


(2.5) 


where  the  primes  have  been  dropped  as  superfluous.  This  relationship  is 

the  same  as  that  developed  for  a propeller  in  a uniform  stream.  Following 

6 

Glauert,  we  now  let 

u)  - 2 a'  n (2.6) 


This  approximation  is  also  inherent  in  the  momentum  theory  in  a uniform  flow 
as  pointed  out  by  von  Mises,°  although  Glauert  is  tacit  on  the  matter. 
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where  a*  is  the  rotational  interference  factor.  Then, 

= 2pn®r^  a'(l-a')  (2.7) 

and  this  expression  for  the  jump  in  pressure  across  the  actuator  disc  must 
be  the  same  as  that  developed  from  the  axial  momentum  theory  as  given  in 
Eq.(l.l2).  Thus,  upon  substituting  Eq.(l.27),  we  find  the  following  relation- 
ship between  the  rotational  and  axial  interference  factors: 


• Lu_ 


a + a 


j*  a rjdr^  = tJPr^a' (I -a') 


(2.8) 


This  relationship  differs  from  the  corresponding  relation  for  uniform  flow 
in  two  respects:  The  terms  U/U  in  the  first  term  on  the  left-hand  side 

CO 

accounts  for  the  parametric  variation  in  the  velocity  across  the  wake, 

d(U/Uo9) 

while  the  term  proportional  to  - - accounts  for  the  shear. 

The  torque  may  be  obtained  by  considering  the  mass  flow  multiplied  by 
the  moment  of  the  rotational  velocity,  thus 


Qa  = Pj  (U+Ujj)w  rd  S 


(2.9) 


Upon  substituting  relationships  already  estab I i shed  or  defined,  there  is 
obtained 


Qg  = hTT  p n J*  (^U  + -yja'r^dr 


(2.10) 


The  apparent  efficiency  is  the  ratio  of  apparent  work  to  the  power 
required  to  generate  it,  i.e.. 


U T 

T)  = 

'a  n 0. 

a 


(2.11) 


We  will  now  determine  the  condition  to  be  satisfied  for  the  propeller 
of  maximum  efficiency,  i.e.,  we  will  minimize  the  power  for  a given  thrust. 
Thus  we  set  the  first  variation  of  flQ  +XU_T,  equal  to  zero  where,  again, 
X is  a Lagrange  multiplier.  For  this  purpose  we  use  the  form  (1.20)  for 
the  apparent  thrust.  Then 

r R Uo,  . R 

n6(ia+^U*®''’a  " ^U+  ^jr^6a'dr  + ^ J r®a'6u^drj- 


+ 2'nXU^p  j (V  + u^)6u^rdr  = 0 


(2.12) 
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We  set  u =2U  a , as  before;  then  by  virtue  of  Eq.(2,8),  we  can  establish 

OD  CD 

the  following  relationship  between  6a'  and  6a 


" r J 6aridri]  = n^  r2(l-2a')6a'  (2.13) 


Upon  substituting  into  Eq.(2.12),  we  obtain 
R 


J r U.=  dr  * 2.)da  - f ^ J da..dr,} 


00  o 


R R 

+ n^U„  J r^a'6adr  + J (V + 2U^a) 6ardr  = 0 . (2.14) 

o o 


Upon  interchanging  the  order  of  integration  in  the  second  term,  we  find 


(U+U  a) 

CD 

l-2a' 


R 

I 

r 


(U+U„a) 

l-2a' 


U = 


i.(iLi 

dr,  \U  / 
1 00 


dri 


+ f2^U„r^a'  + kU^  {V  + 2U^a)]-6ardr  = 0 . (2.15) 

Then,  upon  invoking  the  fundamental  lemma  of  the  calculus  of  varia- 
tions, it  is  found  that  the  condition  for  maximum  apparent  efficiency  is 


CPr^  , 

+ — 5“  a 

uj 


+ X 


2^-^+  2a 

CD  00 


= 0 


(2.16) 


This  equation  is  a generalization  to  shear  flow  of  the  condition  for 

maximum  efficiency  obtained  by  Betz  over  fifty  years  ago.  It  can  be  seen 

that  it  reduces  to  Betz' condi tion  in  the  case  of  uniform  flow.  For  by 

setting  U=Un=U  =constant  , it  becomes 
R “ 


( 1+a) ( l + 2a)  ^ ^ r^  ^ , 
l-2a'  U ^ 

m 


+ \(l+2a)=  0 


which,  upon  dividing  by  (l+2a),  is  identical  to  Betz' cond i t ion  as  quoted  by 
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Glauert  if  we  set  \=-C  . As  in  the  case  of  axial  momentum,  the  presence 
of  the  integral  terms  means  that  the  optimum  distribution  of  the  inter- 
ference factors  depends  not  only  on  local  conditions  but  also  on  conditions 
elsewhere  on  the  actuator  disc,  and  this  dependence  is  caused  by  the  effect 
of  shear. 

In  order  to  determine  the  axial  and  radial  interference  factors  for 
maximum  apparent  efficiency,  it  is  necessary  to  solve  Eqs.(2.l6)  and  (2.8) 
simultaneously.  These  constitute  a pair  of  simultaneous  nonlinear  integral 
equations.  The  pair  can  be  converted  into  a two-point  boundary  value 
problem  by  using  the  definition  (1.6)  in  (2.8)  and  differentiating  (2.16) 
with  respect  to  r . One  condition  to  be  satisfied  may  be  obtained  by 
setting  r=R  in  Eq.(2.l6)  and  another  by  setting  r=0  in  Eq.(2.8).  Tech- 
niques for  solving  two-point  boundary  value  problems  are  readily  available 

q 

(see,  e.g.,  Roberts  and  Shipman^)  and  they  can  be  applied  to  yield  the 
solution.  It  would  appear  to  be  simpler,  however,  to  obtain  solutions  di- 
rectly from  the  integral  equations  themselves.  Since  the  shearing  term  is 
likely  to  be  small,  the  following  scheme  has  been  employed: 

First,  assume  a value  for  X.  (it  can  be  seen  from  the  limit  of  uniform 
velocity  that  \ is  negative).  Next,  solve  the  pair  of  equations  as  simul- 
taneous algebraic  equations  at  each  value  of  r/R  ignoring  the  shear  terms. 
Use  the  values  of  a and  a'  thus  obtained  to  determine  the  shear  terms 
and  solve  the  equations  again  including  the  shear  terms.  Repeat  this  pro- 
cedure until  the  results  converge.  A computer  program  has  been  prepared 
based  on  this  iterative  procedure. 

The  thrust  coefficient  is  to  be  calculated  using  Eq.(1.30),  while  the 
torque  coefficient  is  found  from  Eq.(2.10)  which,  in  coefficient  form,  becomes 


(2.17) 


The  apparent  efficiency  is  then  given  by  Eq.(1.32),  repeated  here  in 
slightly  different  form: 


2TTK., 


(2.18) 
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III.  RESULTS  OF  A SAMPLE  COMPUTATION 

A computer  program  for  determining  the  optimum  distribution  of  the 
axial  and  rotational  interference  factors  has  been  prepared  and  applied  to 
a particular  measured  wake  determined  from  model  tests.  The  ratio  of  stream 
velocity  to  forward  speed  in  the  wake  is  presented  in  Table  1 and  a listing 
of  the  program  in  Fortran  IV  is  given  in  the  Appendix. 

TABLE  I 

MODEL  WAKE  VELOCITY  DISTRIBUTION 


r/R 

CD 

0. 

.225 

.05 

.2142 

.10 

.261 

.15 

.273 

.20 

.307 

.25 

.333 

.30 

.359 

.35 

.382 

.uo 

.413 

M 

.440 

.50 

.473 

.55 

.508 

60 

.541 

.65 

.577 

.70 

.610 

.75 

.640 

.80 

.671 

.85 

.696 

.90 

.720 

.95 

.743 

1 .00 

.768 

TP 

supplied  by  the  sponsoring  agency. 
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Self-powered  model  tests  had  been  carried  out  and  it  was  found  that  this 
wake  was  appropriate  for  the  advance  ratio  J^=.891  and  that  the  thrust 
coefficient  K^=0.17.  In  solving  the  pair  of  equations  (2.8)  and  (2.16), 
it  is  not  possible  to  specify  K.J.,  but  instead  the  Lagrange  multiplier  \ 
must  be  specified;  consequently,  a value  of  \ was  assumed  and  the  result- 
ing value  of  K.J.  computed.  Then  \ was  adjusted  using  a manual  search 
procedure  until  K.|.  became  equal  to  0.17.  It  is,  of  course,  possible  to 
automate  this  search  procedure  but  for  the  limited  number  of  computations 
carried  out  here,  this  would  have  been  unnecessarily  elaborate. 

The  pair  of  equations  (2.8)  and  (2.16)  were  solved  a second  time  with 
the  shear  terms  omitted  and  2U-Uo  replaced  by  U . This  corresponds  to 
considering  U as  a function  of  r only  parametrically,  and  thus  in  this 
case  all  effects  of  shear  are  suppressed.  A comparison  of  the  results  both 
including  and  excluding  the  effects  of  shear  are  shown  in  Figures  1,  2,  and 
3.  The  distribution  of  thrust  loading  was  calculated  using  Eq.(1.34). 

It  can  be  seen  from  Figure  1 that  the  distribution  of  the  rotational 
interference  factor  is  strongly  dependent  on  shear  only  for  the  inboard 
sections  which  do  not  carry  a significant  portion  of  the  load.  On  the  other 
hand,  from  Figure  2 the  optimum  distribution  of  the  axial  interference  factor 
is  strongly  dependent  on  shear  for  the  outboard  sections  and  suggests  that 
a wake-fitted  propeller  designed  without  accounting  for  shear  is  pitched  at 
the  wrong  angle  to  obtain  maximum  efficiency.  Indeed,  the  optimum  load  dis- 
tribution as  presented  in  Figure  3 shows  that  the  tips  should  be  more  heavily 
loaded  when  shear  is  accounted  for  than  when  it  is  not. 

The  measured  model  wake  given  in  Table  1 was  processed"  to  determine 
the  wake  on  the  full-scale  vehicle  and  the  full-scale  wake  data  are  presented 
in  Table  2.  It  was  assumed  that  the  advance  ratio  and  thrust  coefficient  are 
the  same  at  full-scale  as  for  the  model  tests  and  plots  of  the  rotational  in- 
terference factor,  the  axial  interference  factor  and  the  optimum  load  distri- 
bution for  these  conditions  are  presented  in  Figures  4,  5 and  6,  respectively, 
both  with  and  without  shear.  The  relative  effects  of  shear  are  seen  to  be 
similar  to  the  results  for  the  model  wake. 


by  the  sponsoring  agency. 
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TABLE  2 

FULL-SCALE  WAKE  VELOCITY  DISTRIBUTION 


r/R 

u/u 

00 

■ 

0. 

.270 

.05 

.291 

.10 

.314 

.15 

.328 

.20 

.369 

.25 

.394 

.30 

,420 

.35 

.442 

.40 

.473 

.45 

.500 

.50 

.533 

.55 

.567 

.60 

.600 

.65 

.636 

.70 

.669 

.75 

.700 

.80 

.731 

.85 

.756 

.90 

.779 

.95 

.799 

1.00 

.820 

It  is  also  interesting  to  compare 

the  optimum  di 

model  wake  with  the  optimum  distributions  at  full-scale  including  shear  in 
both  cases.  These  comparisons  are  shown  in  Figures  7,  8,  9.  From  Figure  7 
we  see  that  the  rotational  interference  factors  are  virtually  the  same 
between  model  and  full-scale  at  all  radii.  From  Figure  8 we  see  that  the 
axial  interference  factor  has  a small  but  significant  difference  at  all 
radii  and  this  suggests  that  even  when  shear  is  accounted  for,  the  design 
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model  prototype  in  order  to  optimize  its  performance.  This  is  borne  out  in 
Figure  9 where  a comparison  of  the  optimum  load  distribution  for  model  and 
full-scale  indicates  that  the  tips  should  be  more  heavily  loaded  at  the  ex- 
pense of  inboard  stations  in  full-scale  in  comparison  with  model  scale. 

Since  the  momentum  theory  does  not  allow  the  load  to  drop  to  zero  at 
the  tips  as  it  must,  these  conclusions  may  require  some  modification.  Never- 
theless, momentum  theory  in  combination  with  blade  element  theory  has  served 
as  a fairly  reliable  design  tool  in  the  past  and  the  conclusions  drawn  here 
are  very  likely  to  be  qualitatively  correct.  It  is  planned,  as  mentioned  in 
the  Introduction,  to  examine  a more  sophisticated  theory  for  a special  family 
of  wake  profiles  in  which  the  tip  effect  including  shear  is  accounted  for. 

Finally  it  should  be  remarked  that  the  optimum  apparent  efficiency  at 
full-scale  is  smaller  than  at  model  scale  and  that  the  effect  of  including 
shear  is  to  reduce  the  optimum  efficiency  from  that  which  would  be  obtained 
if  shear  were  neglected. 
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APPENDIX 

Given  below  is  a list  of  computer  symbols  and  corresponding  mathematical 
symbols  for  the  program  used  to  solve  Eqs.(2,8)  and  (2,l6)  simultaneously: 


COMPUTER  SYMBOL 


MATHEMATICAL  SYMBOL 


AJ 

initial  value  of 

'^A 

PNI 

initial  value  of 

\ 

NAAJ 

number  of  values 

NNP 

number  of  values 

of  K 

DAJ 

increment  of  J. 

A 

DPN 

increment  of  X 

UR(1) 

values  of  U/U  gi 

(B 

\ ven  i 1 

URT 

value  of  U/U  at  r/R  • 

H 

increment  of  r/R 

URP(I) 

d(U/U  )/d(r/R) 

00 

XI 

a 

X2 

a' 

NOTES  ON  THE  PROGRAM 


Q.TFE 

DET3 

SUBROUTINE  PGl 

BOT 

TOP 


IBM  36O-SSP  integration  subroutine 
usi ng  trapezoid  rule 

IBM  36O-SSP  numerical  differentiation 
subroutine  using  central  differences 

subroutine  to  solve  Eqs.(2.8)  and  (2.16) 
simultaneously  at  given  value  of  r/R  for 
given  values  of  BOT  and  TOP 

integral  term  in  (2.8) 

integral  term  in  (2.16) 
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LISTING  OF  COMPUTER  PROGRAM 

.TY  Pn2F 

DIMENSION  UR(21  )»URP(21  )#X1  C21  >«)C2(21  )»TDP<2I  )«B0T(2I ) 
DIMENSION  FT(21  ),FB(21  ) 

DIMENSION  FKJ(21  )»FK0<21  ) 

DINENSION  CTF<21 ) 

DOUBLE  PRECISION  SHIP 
COMMON  URT»PI.PN»PJ42 
DATA  X1»X2/A2*0.0/ 

DATA  PI/3.1415926/ 

CALL  IFILE<21#’Grj00I’) 

ACCEPT  105»SHIP 
105  FORMATCAR) 

ACCEPT  100.AJ 
ACCEPT  100»PN1 
ACCEPT  101*NAAJ.NNP 
ACCEPT  100#DAJ 
ACCEPT  100#DPN 
101  P0RMAT(21) 

READ(21»100)  <LIRtl).I  = 1.21  ) 

IFT=UR<21) 

100  FORMATCF) 

H=.05 

CALL  DEr3CH.LF.LRP.21.IER) 

IF(IEK.NE.0)  CALL  EXIT 
TYPE  950#SHIP 

950  PORMATf ‘1 •*5X,’0PTIMLM  LOADING  WITH  SHEAR  - 'ABZ/IX# 

1 • JA’RX’LA'RX'KT’RX*K0’RX’Cr*7X'ETA*) 

DO  991  NAJ=1,NAAJ 
PJ4=PI/A.J 
PJ42=PJ4*PJ4 
PN=PN1 

DO  990  ^P=1»N^P 
XI  (21)=. 5 
X2(21 )=0. 

DO  93  L00P=1»10 
DO  30  K=1.21 

FT(K)=<UR<K)*X1 (K))*URP<K)/C1 .-2.*)G(K)) 

XK=K 

FB(K)=X1 (K>*(XK-1 . )*.05 

30  CONTINUE 
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CALL  0TFE<H,FT*T0P.21) 

CALL  0TFE<HfFB»B0T*2I ) 

DO  A0  I=J»21 

40  TOP(I)=T0P<2I >-TOP<I) 

Y1=X1 (21 ) 

Y2=X2(2I) 

DO  50  J=l#21 

I=22-J 

XI=I 

RI=H*<XI-1 .) 

(^L  PGl  (YI  >Y2»UR<  I )»URP(I  )»TOP<  I )/BOT<I  >»RI  > 

Xl(I)=Yl 

' X2(n=V2 

' 50  cofsrriNUE 

90  CONTINUE 

DO  R0  I=l»21 
xr=T 

R=(XI-1.)*H 

FK.  I(  T )=(  (2.*UR(  I >-LRT)+Xl  ( I ) )*X1  < I )*R 
FKO( I )=(UR( I )+Xl ( I ) )*X2 ( I )*R*R*R 
IF  (R.EO.0*)  (DTO  81 

CTF( T )=8.*R* (UR( T )*X1 ( I )*X1 < I )*X1 < I ) >-8«*URP( I )*BOT< I ) 
GOTO  80 

81  CTF<I>=0. 

80  CONTINUE 

200  FORNWT<2X»F4.2j4E16.4) 

CALL  0SF(H»FKJ>T0P*21 ) 

CALL  0SF<H*FKQ»B0T»21> 

CT=R.*T0P(21 ) 

Fvl=PI*AJ*AJ*TOP(21  ) 

F0=.5*PI*PT*A>J*B0T(21  ) 

EH=.5*AJ*FJ/F0/PI 
TFCNP.GTT.l)  GOTO  95 
TYPE  900#AJ»PN»FJ<FO»CT#EH 
GOTO  990 

95  TYPE  905#PN*FJ#FO#CT»EH 

900  FOR'^T(2X»6F10.4) 

905  FORMAT<12X»5F10.4) 

99(^  PN=PN-DPN 

991,  A.I=AJ*I>AJ 

END 
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SUBROUTINE  PGl ()U#X2»UR»URP#TCP*B0T»RI ) 
COMMON  URT>PI»PN»PJ42 
DO  10  1=1 >50 

F1=(UR*X1  )*(LR+2.*X1  )-T0P*<l.--2.*X2) 
Fl=Fl*PJ4e*RI*RI*)C*<l .-2.*)C) 
F1=F1*PN*(1 .-2.*X2)*(2.+UR-URT+2.*X1 ) 
IF<RI.EO.0.)  GOTO  30 

F2=UR*X1  ♦Xl  ♦Xl -URP>fBOT/RI 

GOTO  40 

30  F2=UR*X1+X1*X1 

40  CONTINUE 

F2=F2-PJ4e*RI*RI*X2*(l .-X2) 

FSO=Fl *F1  ♦F2*F2 
IF(FS0.LT.1E-12)  RETURN 
F1X1=3.*UR^4.*X1+2.*PN*<1.-2.*X2) 

FI  X2=2  . *TOP+P,  •^♦R I *R  I * ( 1 . -4.  *X2  ) 

FI  >C=F1  X2-2.*PN*<2.*UR-URT-»2.*X1  ) 
F2XI=UR+2.*X1 

F2)C=-PJ42*RI*RI*<1 .-2.*X2) 

XI=F1 XI *F2X2-F1 X2*F2X1 
X1=X1*(F2*F1X2-F1*F2X2)/X.) 
je=X2+(Fl*F2Xl-F2»FlXl )/XJ 
10  CONTINUE 

RETURN 
END 
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FIGURE  3.  THRUST  LOADING  DISTRIBUTION  WITH  AND  WITHOUT  SHEAR 
MODEL  WAKE 


FIGURE  7.  ROTATIONAL  INTERFERENCE  FACTOR  DISTRIBUTION  WITH  SHEAR 
IN  MODEL  AND  FULL-SCALE 


IN  MODEL  AND  FULL-SCALE 
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FIGURE  9.  THRi 
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